Influence of strong microwave radiation on static properties of a small Josephson 

junction 
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Preliminary results are presented concerning static properties of a small Josephson junction under 
the influence of strong microwave radiation. We discuss the correspondence between a Brownian 
particle moving in a periodic potential and superconducting phase difference in a small Josephson 
junction. Next, we describe an experimental method of determining the amplitude of microwave 
current flowing across the junction. Typical examples of static characteristics of the junction are 
presented, including its dynamical resistance as a function of microwave power. We discuss also 
the influence of an external magnetic field on the junction dynamics and show that in this case the 
one-dimensional Stewart-McCumber model becomes insufficient. 

PACS numbers: 74.50. +r, 05.45.-a 



I. INTRODUCTION 

Recently, the noise-assisted transport of Brownian 
particles has attracted considerable attention, which is 
still growing in view of possible applications in various 
branches of physics and chemistryjir— Unfortunately, ex- 
perimental investigation of transport phenomena on a 
molecular scale is extremely difficult, thus apart from 
theoretical studies there is a need for adequate physical 
systems, modeling various transport phenomena in a way 
which could be accessible experimentally. 

As shown in Refs.|j|a, there is a rigorous mathematical 
equivalence between the Newton-Langevin (NL) equa- 
tion, describing a Brownian particle in a periodic po- 
tential, and the Stewart-McCumber (SM) model deal- 
ing with the superconducting phase dynamics in a small 
Josephson junction in the presence of thermal fluctua- 
tions. In particular, a constant external force acting on 
a Brownian particle corresponds to a constant current 
biasing the Josephson junction, while a time-averaged 
velocity of the particle is equivalent to a constant volt- 
age induced on the junction electrodes. In other words, 
the relevant parameters of a Brownian motor can be re- 
lated to easily measurable macroscopic quantities such 
as a constant current and constant voltage, leading di- 
rectly to the current-voltage [I-V) characteristic of the 
Josephson junction. 

The aim of the present paper is twofold. First, we pro- 
pose a uniform method making possible to determine the 
absolute value of microwave current amplitude between 
junction electrodes, independently of the junction geom- 
etry and its coupling to the experimental setup. Such a 
uniform approach covers a wide range of junction param- 
eters, including two limiting cases of a "current source" 
and a "voltage source" model. Second, we investigate 
the influence of a strong microwave signal on the static 
properties of a small Josephson junction, with particular 
attention paid to its dynamic (differential) resistance. In 
this context we present preliminary results and explore 
the possibility of detecting a negative resistance (both 
absolute and differential) on the I-V characteristic. 



The paper is organized as follows. In Sec. II we discuss 
the equivalence between the NL and SM models. Next, 
we recall the well-known relation between the critical cur- 
rent I{) and the normal resistance i?„ of a real Josephson 
junction, and discuss the possibility of making the junc- 
tion parameters compatible with theoretical estimates. 
Experimental details, including coupling of the junction 
to a microwave setup, are given in Sec. III. Sec. IV 
deals with an experimental evaluation of rf (microwave) 
current flowing across the junction. In particular, we 
propose an analytical approximation making possible to 
determine the rf current amplitude for a wide range of 
junction parameters. A few examples of the current- 
voltage [I-V) characteristics are presented in Sec. V with 
a special attention paid to the dynamic resistance as a 
function of strong microwave signal. We discuss also the 
influence of an external magnetic field on the junction 
behavior and show that the superconducting phase dif- 
ference (/) becomes spatially modulated, making the SM 
model insufficient for a proper description of the junction 
dynamics. Sec. VI contains concluding remarks, we in- 
dicate also possible generalization to a multidimensional 
model of a Josephson junction immersed in an external 
magnetic field. 



II. MODELING OF CHAOTIC PHENOMENA IN 
A JOSEPHSON JUNCTION 

The chaotic dynamics of a classical Brownian particle 
in the presence of thermal noise can be described by the 
inertial NL equation)^ 



mi + -ii + V'{x)=F + Acosm+y/2^^k^i{t) (1) 

where a dot and prime denote differentiation with respect 
to t and X, respectively. The particle mass is denoted by 
TO, 7 is the friction coefficient, $7 — the angular frequency, 
F and A denote constant and alternate external forces, 
respectively. Thermal fiuctuations are described by zero- 
mean Gaussian white noise ^(i) with the autocorrelation 
function {^{t)£,{u)) — 5{t~u), fc^ denotes the Boltzmann 
constant and T — temperature. 



If a spatially periodic potential V{x) is cosinusoidal, 
then it can be easily verified that Eq. (1) is formally 
equivalent to the SM model describing the supercon- 
ducting phase dynamics of a one-dimensional small-area 
Josephson junction*^'^ 



Ccl) 



Rn 



(jy + h sin((/)) 



= Id + Ia COs{nt) + 



2kBT 

Rn 



m 



(2) 



where cj) denotes the superconducting phase difference be- 
tween electrodes and the junction is characterized by its 
critical current /q, capacitance C and normal-state resis- 
tance Rn ■ The amplitudes of the dc and ac (microwave) 
currents applied to the junction are denoted by Id and 
la, respectively. 

A dimensionless form of Eq. (2) can be written as 



di'2 



-f CT-^ + sin((/«) =iQ + ii cos(f^ii') + V2aDr{t') (3) 



dt 



where the time has been rescaled to t' ~ i/io, tq = 
l/u!p — y^hC/2eIo- Consequently, fli — flro, a ~ 
T^/RnC, and the noise intensity D = 2ekBT /TiIq. The 
dimensionless dc and ac amplitudes are given by iq — 
Id/h and ii = la/h, respectively. 

Comparison of Eq. (3) with Eq. (1) shows clearly that 
(after appropriate scaling of Eq. (1)) both models are 
mathematically equivalent, although the physical mean- 
ing of relevant variables is distinctly different. In partic- 
ular the superconducting phase difference (p in Eq. (2) 
corresponds to the spatial coordinate x in Eq. (1). Con- 
sequently, the average value of the time-derivative {4>t), 
proportional to the constant voltage across the junction, 
is analogous to the averaged velocity of a Brownian par- 
ticle. Similarly, external current amplitudes Id and la 
correspond to external forces F and A, respectively. 

The range of dimensionless parameters in which we 
can expect either an absolute or a differential negative 
resistance, is given in Ref. |5|. Thus, using the scaling 
of Eq. (3) and coming back to Eq. (2) it is possible to 
recover absolute values of the junction parameters and 
the external driving. 

For example, the noise intensity D — 10^"^ implies Iq ~ 
176 ^A at r = 4.2 K. Moreover, assuming fi = 150 x 
10^, i.e. / ~ 23.87 GHz we obtain the following sets of 
junction parameters 



(a) Rn = 2.6 f7, C = 13.8 pF, Iq = 176 ^A 



(b) Rn = 1.4 f7, C = 27.9 pF, Iq = 176 ^A 



(4) 



The external rf current amplitude may vary, ranging 
from la — 120 fiA for the set (a) up to la — 500 ^A for 
the set (b),^ 

The main difficulty with practical application of the 
above parameters follows from the fact that for a real 



Josephson junction the normal resistance i?„ and the crit- 
ical current Iq are not independent, but their product is 
given by:^ 



laRn — 



2e \2kBTJ 



(5) 



where A(r) denotes the superconducting gap energy. 

For niobium at T = 4.2 K, we have A(T) ~ 1.4 meV, 
thus the upper limit for /oi?n can be estimated according 
to Eq. (5) as I^Rn — 2 mV. In practice, the value of 
IqRu for a real junction is lower due to inevitable defects, 
nevertheless it is still a few times greater than the value 
oi IoR„ resulting from the set of parameters (4) required 
by the theory. 

In order to make real junction parameters compatible 
with theoretical predictions one can use one of two ex- 
perimentally accessible methods, reducing either i?„ or 
Iq. In the first case, the effective value of normal resis- 
tance can be decreased by shunting the junction with an 
additional resistor connected in parallel to the junction 
electrodes^ In this work, however, we use an alterative 
approach to reduce the critical current Iq by applying an 
external constant magnetic field in the junction plane. 

The following sets of junction parameters have been 
chosen for further studies: 



(a) Rn = 2.96 n, C ^ 14.0 pF, Iq = 460 ^lA 



(b) Rn = 1.45 n, C ^ 28.0 pF, Iq = 920 /lA 



(6) 



Comparing the above parameters with theoretical pre- 
dictions (4) one can see that the values of i?„ and C 
correspond approximately to those specified previously, 
while the nominal critical currents Iq are a few times 
greater than Iq ~ 176 /lA required by the theory. As 
mentioned above, the critical current can be adjusted by 
applying an external magnetic field to the junction. In 
this manner, we are able to obtain the relevant junction 
parameters in agreement with those predicted theoreti- 
cally. 

The amplitude /„ of the rf current flowing across the 
junction is obviously proportional to vP, where P de- 
notes the incident microwave power. Unfortunately, la 
depends in a rather complicated way on the geometry 
of the sample holder and its coupling to the microwave 
source. Therefore, the absolute value of la has been de- 
termined experimentally by using the well-known rela- 
tions between the effective constant current at zero volt- 
age and the rf current amplitude la (see Section IV for 
details). 



III. EXPERIMENTAL 

The Nb-A10a;-Nb Josephson junctions in overlap ge- 
ometry have been designed and manufactured in the In- 
stitute of Radio Engineering and Electronics in Moscow. 



The sample holder was designed as a section of a ridge 
waveguide of reduced impedance, which was used to im- 
prove matching of the junction to the standard K-hand 
(WR 42) waveguide and obtain a relatively large ampli- 
tude la of the rf current. Additionally, la could be con- 
trolled by a variable short located just behind the sample 
holder. Two superconducting coils of nearly rectangular 
cross-section were attached to the side walls of the sam- 
ple holder, and could be used to apply a dc magnetic field 
of up to ~ 40 Gs directed parallel to the junction plane. 

The measurement system consisted of a _?C-band 
klystron, followed by a ferrite isolator, ferrite modula- 
tor and calibrated attenuator. The microwave frequency 
could be controlled in the range 23 - 25 GHz and the 
maximum output power level was estimated as 100 mW. 

The current-voltage (I-V) characteristics were mea- 
sured by a standard 4-point technique. The junction was 
biased by a digital current source while the voltage was 
measured using a digital voltmeter and recorded by a 
computer. All the measurements were carried out at the 
temperature T = 4.2 K. 




FIG. 1: (Color online) The relative rf current amplitude ii at 
Imax = plotted as a function of the normalized frequency ^. 
Full squares denote numerical results, the dashed line shows 
an asymptotic linear dependence, and the solid line represents 
an analytical approximation (10). 



IV. EVALUATION OF THE MICROWAVE 
CURRENT AMPLITUDE 

According to the "current source" model^ we consider 
the following equation: 



-— + sin((/)) = iQ 
dr 



Zisin(^T), 



(7) 



where io — Id/Io, h — la/h as before, <^ — Sl/wc and 
Wc = 2eRnIo/^- As compared to Eq. (3), we use a 
slightly different time scaling r = flct and neglect both 
the second time-derivative and the noise term. 

Eq. (7) has been solved numerically in several 
papersii2r2i and it has been shown that the maximum 
constant current Imax at zero voltage decays with in- 
creasing rf amplitude ii. However, detailed behavior of 
Imax as a function of zi depends strongly on the fre- 
quency parameter ^. For ^ <C 1, we have a "current 
decreases linearly and Imax = is 
On the other hand, for ^ ^ 1 the 
'voltage source" with a Bessel-like 



source" model, Ima-. 
attained for ii ~ 1. 
model tends to the 
behavior: 



c/Io — Jo 



2eR„Ia 

nn 



(8) 



where Jq denotes the Bessel function of order and the 
argument can be transformed to a simpler form by using 
dimensionless quantities: 



2eRnIa 
nil 









(9) 



Hence, asymptotically Imax = is obtained at ii — 
Ko,iCi where Ko,i = 2.405 denotes the first zero of the 
Bessel function Jq. Fig. 1 shows the results of numerical 



calculations of «i Ir — n as a function of £, taken from 
Refs. Pzllgl-llli The dashed line denotes an asymptotic lin- 
ear dependence ii = ko,iC and the solid line shows an an- 
alytical approximation suggested here to cover the whole 
range of £,: 



«i = Y^i + («o,i6' 



(10) 



It is clear that the above approximation exhibits a 
proper behavior both for ^ — j' and ^ 3> 1. Moreover, 
as shown in Fig. 1, the analytical approximation (10) 
is in good agreement with numerical data, thus it can 
be used to evaluate the rf current amplitude ii also for 
intermediate values of ^. 

As an example. Fig. 2 shows an experimental depen- 
dence of Imax on the relative amplitude of microwave 
signal (P/Pq)^/^, where P/Po denotes the relative mi- 
crowave power. 

The junction parameters are slightly different from 
those given in Eq. (6) and are equal to P„ = 2 il, 
Iq = 680 /iA, hence ^ ~ 0.037 for the microwave fre- 
quency / ~ 24 GHz, and according to the analytical 
approximation (10) we find ii|/,„^^=o — 1.004. It follows 
from Fig. 2 that Imax = for"(P/Po)i/2 ^ o.45 what 
means that the absolute amplitude of the rf current at 
this point is equal to la — Iq = 680 fiA. Thus, the max- 
imum available current amplitude can be evaluated as 
jijnax) _ j^^Q^^g ^ 2.5 niA and is far above the range 
of la predicted by theoretical considerations (see Section 
II). Strictly speaking, the above analysis deals with the 
RS J model, since the second time-derivative has been ne- 
glected in Eq. (7). The influence of the junction capaci- 
tance on the rf current amplitude is not fully understood, 
however as follows from Ref. |Tl|, the estimates taking into 
account a nonzero capacitance yield even larger values of 







• 








1 1 1 1 1 1 1 










• 












6UU 


- 


• 


-• 






- 












~« 




■ 




500 


- 








• v 


- 


^-^ 
















< 


400 


- 








^i^v 


- 


3. 
































S 


300 


_ 








'■• 


_ 


^ 












"•. 


■ 




200 


- 










- 




100 

n 










1 1 1 1 1 ~ 1 


■ 



200 



0.0 0.1 0.2 0.3 0.4 0.5 

FIG. 2: (Color online) Experimental dependence of the max- 
imum constant current Imax (full circles) on the relative am- 
plitude of the microwave signal [P/PoY ■ The dashed line 
shows a linear fit to the experimental data. 
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FIG. 3: (Color online) Family of I-V characteristics plotted 
as a function of incident microwave power level. 



To conclude this section, we were able to estimate the 
rf current amplitude for a wide range of dimensionless 
parameter ^, which depends on the junction parameters 
and the rf frequency. In particular, it was shown that 
the rf amplitudes required by the theory could be easily 
obtained in our experimental setup. 



V. I-V MEASUREMENTS — RESULTS AND 
DISCUSSION 

In order to study the influence of microwave radiation 
on the static properties of a Josephson junction, we have 
performed a series of detailed I-V measurements for var- 
ious levels of microwave power. Twelve sets of junctions 
have been measured with the parameters close to those 
given in Section II. 
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FIG. 4: (Color online) Enlarged part of Fig. 3 plotted for 
smaller values of current and voltage. 



The results obtained so far exhibit a rich variety of 
different phenomena which are still under investigation. 
Here we present only a few examples, dealing mainly with 
the dynamical resistance of the junction in the presence 
of strong microwave radiation. 

A typical example of the I-V characteristic measured 
for the sample (a) is presented in Fig. 3, while Fig. 4 
shows an enlarged part of the same diagram plotted for 
smaller current and voltage values. According to Section 
IV, the power level P/Pq = — 20 dB corresponds to the 
absolute rf amplitude la — 150 fiA and consequently —10 
dB denotes la — 480 /iA. For the set of parameters given 
by Eq. (4), we have i7i < 1 and ili > a. As a conse- 
quence, the Shapiro steps are not observed in this region, 
in agreement with Refs. ISlM 

It is clear that the I-V characteristics are influenced 
strongly by the microwave radiation. First, we observe 
a step in the I-V curve, which corresponds to a local 
decrease of the dynamical resistance Rd = dV/dl from 
an asymptotic value Rd — 5.3 il, to Rd — 1.1 ft — 1.5 SI. 
The voltage position of this step, denoted by Vs , depends 
linearly on the rf amplitude as shown in Fig. 5. 

The second effect of the microwave radiation is its in- 
fluence on the dynamical resistance Rd at zero voltage 
(see Fig. 6). This time we observe rather large and ir- 
regular changes of Rd ranging from a very small value 
Rd < 0.1 n for {P/PoY/^ = 0.177 up to Rd ^ 8.4 ^ 
for (P/Po)^/^ = 0.224. It is clear that the depen- 
dence of Rd on the rf radiation is non-monotonic and for 
(P/PoY^^ = 0.177 we observe nearly zero value of the 
dynamical resistance. Unfortunately, contrary to theo- 
retical predictions, so far we were not able to observe a 
negative dynamical resistance of the Josephson junction. 

It seems that the main reason for a visible discrepancy 
between theoretical and experimental results follows from 
the method used to reduce the critical current of the junc- 
tion. A more detailed analysis shows^^ that applying an 
external magnetic field in the junction plane gives rise 
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FIG. 5: (Color online) Step position Vs as a function of the 
relative rf amplitude. Experimental data are plotted as full 
circles and the dashed line shows a linear fit. 
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FIG. 6: (Color online) Dynamical resistance at zero voltage 
as a function of the relative rf amplitude. 



to a spatial modulation of the current density. In other 
words, the external magnetic field makes the phase (^ to 
become a function of spatial coordinates, in spite of the 
fact that we deal with small junctions whose physical di- 
mensions are much smaller than the Josephson penetra- 
tion depth A J. As a consequence, the Stewart-McCumber 
model becomes inadequate and should be generalized to 
a multidimensional case, taking the form of a partial dif- 
ferential (sine-Gordon like) equation: 



-aV^(/)-h- 



dt 



-sin(0) = io-hiiCos(fit)-hV2aZ?r(i) 

(11) 



where V^(/) = d?(j}/dx'^ + d^(j)/dy'^ denotes the two- 
dimensional Laplacian, a is a scaling parameter and 
4>{x, y, t) is a function of time and two spatial coordi- 
nates. 

Unfortunately, such a model has no counterpart in the 
domain of chaotic dynamical systems, hence the theo- 
retical predictions resulting from the Newton-Langevin 
equation cannot be directly applied to a Josephson junc- 
tion immersed in an external magnetic field. 



VI. SUMMARY AND CONCLUSIONS 

In this work we investigated the infiuence of high- 
power microwave radiation on the static characteristics 
of a small Josephson junction. In particular, the depen- 
dence of the critical current on the microwave power level 
was studied, allowing to determine the absolute value of 
the rf (microwave) current flowing across the junction. 
It was shown that for the junction parameters as given 
in Section II, a "current source" model can be applied, 
leading to a simple relation la — Iq at the point where 
the maximum constant current Imax attains zero. 

In the next step, a dynamical resistance of the junction 
was investigated as a function of the rf current amplitude. 
In particular, it was found that the dynamical resistance 
at zero voltage is strongly dependent on the microwave 
power level and may vary in a wide range, however neg- 
ative values of the resistance were not observed. 

As far as the correspondence between Newton- 
Langevin and Stewart-McCumber models is concerned, 
we found that the junction parameters required by the 
theory were incompatible with those of really existing 
Josephson junctions due to a natural constraint on the 
product of R„ and Iq. As pointed out in Section V, an 
external magnetic field applied in the junction plane may 
reduce the critical current, but also makes the supercon- 
ducting current density to be spatially modulated. In 
other words, the analogy between the NL and SM mod- 
els breaks down, since the dynamics of a Josephson junc- 
tion immersed in an external magnetic field should be 
described by a more general multidimensional equation. 
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